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Die Methode der numerischen Integration findet auf dem Gebiet der Astronomie ihre fruchtbarste 
Anwendung in der Behandlung des Dreikörperproblems und sie ist in dem Falle, als die Voraussetzungen, 
dasselbe als Störungsproblem zu behandeln, das heißt: die Keppler’sche Bewegung als erste Näherung 
ansehen zu können, nicht gegeben sind, der einzige Weg, der in einem konkreten Falle die Bestimmung 
der Bewegungen in dem betreffenden Maßensystem ermöglicht. 

Für das Zweikörperproblem mit dem relativ einfachen funktionellen Zusammenhang der Koordinaten 
und der Zeit scheint von vornherein kein Anlaß vorzuliegen, diese Methode in Anspruch zu nehmen, und 
das um so weniger, als sie gerade bei der Ermittlung von Größen nullter Ordnung viel von ihrer 
Bequemlichkeit verliert. Nun kann man aber zeigen, daß sich bei der numerischen Integration der Glei¬ 
chungen dieses Problems ein ähnliches »direktes« Verfahren angeben läßt, wie es Th. v. Oppolzer für die 
Encke’sche Methode der Ermittlung der speziellen Störungen in den rechtwinkligen Koordinaten getan 
hat. Dasselbe kommt bekanntlich darauf hinaus, daß man durch Verbindung der dynamischen Gleichung 
mit dem Ausdruck, den die numerische Integration für das Doppelintegral gibt, zugleich den Funktions¬ 
wert und den Wert des Doppelintegrals für dasselbe Argument bestimmt und erst nötig hat, die zweiten 
Differenzen zu extrapolieren, was ja immer mit ausreichender Sicherheit geschehen kann, so daß die 
Bezeichnung »direktes Verfahren« praktisch völlig gerechtfertigt erscheint. Ist aber ein derartiges Ver¬ 
fahren beim Zweikörperproblem möglich, so ist dadurch ein Mittel gegeben, die Koordinaten auf eine 
sehr rasche und einfache Weise zu finden, das in dieser Hinsicht die gebräuchlichen Methoden unter 
Umständen wesentlich übertrifft. 
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Es scheint mir daher diese Anwendung für die Ephemeridenrechnung von praktischer Bedeutung 
zu sein, da damit gerade der umständlichere Teil derselben, wie die Auflösung der Kepplerschen 
Gleichung, zu Gunsten eines weit einfacheren Verfahrens entfällt, ein Umstand, der bei stark exzen¬ 
trischen, aber noch nicht parabelnahen Bahnen besonders zum Ausdruck kommt. 

Die hier notwendigen Formeln der numerischen Integration sind darin von den gebräuchlichen 
verschieden, daß die Ausgangswerte des einfachen und Doppelintegrals von Null verschieden sind. Es 
sollen daher zunächst die diesbezüglichen Formeln entwickelt werden. 


I. 


Die Bessel’sche Interpolationsformel lautet in der bekannten Gauß-Encke’schen Bezeichnungsweise 
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wo Ei+ 1 / 2 der Wert des einfachen Integrals für in = 0, das heißt für das Argument a + [/ + —) co, ist. 


Durch nochmalige Integration erhält man 
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wo + i/ 3 wieder der Wert des iterierten Integrals für das Argument a + ^ ] <o ist. Für die Grenzen 

1 1 

m — -bis 4-ist dann 

2 2 
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u. s. w. ist. 

Summiert man über sämtliche i von 0 bis i— 1, so wird 
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Ist weiter £ 0 der Wert des einfachen Integrals für l — a, so ist 
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daher nach den bekannten Entwicklungen der numerischen Integration 
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Bestimmt man das willkürliche Anfangsglied l f\ a -—j aus 

V 2 1 
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so reduziert sich der Ausdruck für— Ei.+ q s auf die erste Zeile und die Summe 


co 


- [Eij 2 4- Ezj z 4- .... 4- Ei— i/J 


wird 


= H/(a + i (ü) - u /(a) + Pl [f(a + i co) -f(a)] + p 2 [/“ (a + i co) -/» (a)] + 
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Demnach ist 


^ JJ f (0 dl 2 — ll f(ci + i w) — u f (a) + (p l + cj [f(a 4- i w) —/(«)] 4- 


+ (ft + <*) [/" (a + *«0 -/ n («)] + 


Bezeichnet man mit J t und J 0 den Wert des Doppelintegrals für das Argument a 4 - zoo, beziehungs¬ 
weise a , setzt weiter 
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so wird 
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CO 2 
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CO 2 

ein Ausdruck, der sich wieder auf die erste Zeile reduziert, wenn der willkürliche Ausgangswert der 
zweiten Summenreihe 

"/(«) = 4- ßi/(«)- ö 2 / n («)-•••■ 

OO 2 


gesetzt wird. 

Zur numerischen Bestimmung des Doppelintegrals hat man daher folgendes Rechnungsschema: 

Sind E 0 und J 0 die Werte des einfachen, beziehungsweise iterierten Integrales für das Ausgangs¬ 
argument a und bestimmt man die Anfangswerte der ersten und zweiten Summenreihe aus 
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so ist der Wert des Doppelintegrales J,- für das Argument a 4- ita gegeben durch die Formel 
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il. 


Die Methode der numerischen Integration kann zunächst auf eine Gleichung angewendet werden, 
die den Radiusvektor als Funktion der Zeit definiert. 

Führt man in die Gleichungen des Zweikörperproblems 


d 2 x 79 x - 
— + k 2 — = 0 , 
dt 2 r 3 

die Polarkoordinaten r und v ein, so wird daraus 


dt 2 r 3 



d 2 r (d vY 



d r 

d v 

d 2 v1 

COS V 

— — r — 

— sin v 

2 

— . 

— 

+ r — 


dt 2 \dt) 



d t 

dt 

dt 1 


\d 2 r (dvY\ 


r 

dr 

d v 

d 2 v 1 

sin v 

- -r — 

+ COS V 

2 

— , 

,- 

+ r — 


dt 2 \dt) _ 



d t 

d t 
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k 2 
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woraus einerseits 


d 2 r fdv\ 2 k 2 ^ 
- rt — + - = 0 


dt 2 


d t 


dv 

andrerseits der Flächensatz folgt. Eliminiert man mit Hilfe des letzteren—, so erhält man 

d t 


d 2 r k 2 p k 2 

-1- 

dt 2 r 3 r 2 


also eine Gleichung in r allein. 

Denkt man sich nun für eine Reihe von äquidistanten Argumenten t die Werte der Funktion 

d 2 r 

~~dt 2 

gegeben und die Differenzen- und Summenreihen in entsprechender Weise gebildet, so ergibt sich der 
Wert des Doppelintegrals r aus 
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wo mit o) wieder das konstante Argumentintervall bezeichnet sein soll. 


d 2 r 


Im Verein mit der obigen dynamischen Gleichung läßt sich daraus sowohl — als auch r bestimmen, 
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wenn die Größen IX /,/ n ,.. . bekannt sind. Angenommen, es sei —und r für das Argument t 0 + (i —1) co 

dt 2 

bekannt, so folgt daraus der Wert ll f ( t 0 + i w) und damit auch die einzige Größe, die zur Auswertung von 
d 2 r 

— und r für das nächste Argument t 0 + z*to notwendig ist. Bei dem relativ engen Intervall einer Epheme- 
d t~ 

d 2 r 

ridenrechnung werden die zweiten Differenzen von — mit Rücksicht auf den kleinen Koeffizienten wohl 

dt 2 

immer zu vernachlässigen oder in besonderen Fällen mit ausreichender Sicherheit zu extrapolieren sein. 

d 2 r 

Die Elimination von —- aus den beiden Gleichungen ergibt unter diesen Umständen 
dt 2 
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r — to 2 u f -+- 


m 2 k 2 p—r 
12 ' r 3 


zur Bestimmung des Radiusvektors. Die Auflösung kann im allgemeinen auf sehr einfache Weise 

CO2 W Y 

bewerkstelligt werden - . -bedeutet ja für gewöhnlich einen sehr kleinen Betrag. Für Planeten 

12 r 3 

ist das unmittelbar ersichtlich, gilt aber auch für Kometenbahnen, wenn die Periheldistanz q nicht aus- 
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nahmsweise klein ist. - -hat im Perihel den Wert —---, nimmt zunächst ab, geht für r—p durch 


Null und erreicht bei r 
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— p den größten negativen Wert, dessen absoluter Betrag— . — ist. Letzterer 
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ist aber nur für s < 0‘ 11843. . . größer als im Perihel, so daß für Kometen —-— der größtmögliche 
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1 p —Y 

Betrag bei gegebener Bahn und - der größte Betrag überhaupt ist, den~ -bei einer vorgelegten Perihel- 

q 1 r 3 

ö) 2 k 2 

distanz q annehmen kann. Nun ist — . eine Größe, die bei mäßigen Intervallen und nicht extrem 

12 q 2 

kleinen Periheldistanzen stets klein sein wird. Für q = 0'l, eine Sonnennähe, die nur bei etwa—der bis- 

30 


her bestimmten Kometenbahnen beobachtet wurde, ist bei viertägigem Intervall 

0t)2 P _ Y 

man demnach in- . -—— für r als erste Näherung co 21 ^ so wird 
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= 0*04. Setzt 
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unter gewöhnlichen Umständen von ausreichender Genauigkeit sein und sich nur selten die Notvvendig- 

0t)2&2 p _ y 

keit ergeben, diesen Wert neuerdings in-. - - einzuführen. Es tritt hier noch der günstigeUmstand 

12 r s 

dazu, daß diese Annäherung um so besser ist, je kleiner die Exzentrizität ist. Für parabolische oder 
parabelnahe Bahnen hat man aber genügend rasche Methoden der Koordinatenbestimmung, als daß ein 
Anlaß vorliegen würde, andere rechnerische Grundlagen aufzusuchen; es wird also eine der Einheit nahe 
Exzentrizität hier ohnehin nicht in Betracht kommen. Die vorliegende Methode wird überhaupt in jenen 
Fällen ihre besondere praktische Verwertung finden, in denen weder die Auflösung der KeppleFschen 
Gleichung noch der Anschluß an die Parabel ein bequemes Rechnungsverfahren bedeutet. 


Man bestimmt weiter 

— = k 2 . v ~ r 
dt 2 ‘ r 3 

oder mit genügender Annäherung 

d 2 r _ jp — to 211 / 

"dt 2 w 611 / 3 

also den Wert der Funktion f für das Argument t 0 + im, woraus n f[a + (i + l)o>] erhalten wird und 
damit r für das folgende Argument a + (i + 1) w. 

Man kann übrigens von vornherein die Annäherung um einen Grad weiter treiben. Da das ent¬ 
sprechende Verfahren eine kaum nennenswerte Mehrarbeit bedeutet, ja sogar unter allen Umständen 
eine zweite Durchrechnung überflüssig macht, außerdem aber den Grad der Annäherung und damit das 
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numerische Ausreichen dieser Methode direkt beurteilen läßt, so wird es sich empfehlen, diese weitere 
Annäherung gleich von Anfang an einzuführen. 


Es ist in 


d ^L — k 2Pzzl 


dt 2 


für r zunächst co 2li / substituiert worden. Setzt man dafür den verbesserten Wert 
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daher 
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Was die Größe v anbelangt, so genügt für diese eine ganz rohe Annäherung, die man immer sofort 
hersteilen kann. 
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In der Umgebung des Perihels, die man mit besonderem Vorteil als Ausgangsstelle der Ephemeriden- 
rechnung nach dieser Methode wählen wird, ist 
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Ist die Rechnung aber bereits im Gange, so ist für die Ermittlung von v mit genügender Annäherung 
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Mit diesem so verbesserten — wird man schließlich 

dt 2 


r — (ü 2l] f + 


co 

12 


<2 2 r 


rechnen und damit einen Wert erhalten, der auch für eine strenge Ephemeridenrechnung völlig ausreicht, 
wenn nicht ganz außergewöhnliche Sonnennähen stattfinden. 

Es erübrigt nun noch, die zur Bildung der beiden Summenreihen notwendigen Ausgangswerte r 0 , 
ld 2 r\ 

und — zu bestimmen. Bezieht sich die Ephemeride auf ein Bahnstück, das das Perihel nicht 

\ d i 2 ) 0 

enthält, auch nicht in der Nähe desselben gelegen ist, so hat man für die gewählte Ausgangsstelle r 0 , 
beziehungsweise v Q zu ermitteln und erhält dann 
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ln 


Da dieser Fall wohl nur bei Oppositionsephemeriden von Planeten eintreten wird, bei denen wegen 
der Kleinheit der Exzentrizität die gewöhnliche Methode keine Unbequemlichkeit in sich schließt, so wird 
man sich die Ausgangswerte nach dieser leicht verschaffen können, wenn man das hier gegebene Ver¬ 
fahren auch bei kleinen Exzentrizitäten anwenden will, wo es allerdings kaum einen wesentlichen Vor¬ 
teil gewährt. Enthält aber der Bogen, auf den sich die Ephemeride bezieht, das Perihel, was bei Kometen 
wohl stets der Fall ist, so können diese Ausgangswerte in sehr einfacher Weise direkt aus den Elementen 
gewonnen werden. 


Im Perihel selbst ist 


r P = q, 




sk 2 {d 3 r\ _ fd 4 r\ 
q 2 ’ [ dt 3 ) P ~ ’ [dttjp 


s ( 3s + l) /fe 4 

r 


Wählt man nun als Ausgangsepoche etwa die dem Periheldurchgang benachbarte mittlere Mitter¬ 
nacht des Meridians, auf den sich die Ephemeride bezieht, und ist A t der Abstand dieser Epoche vom 
Perihel in mittlerer Zeit, so ist mit genügender Annäherung 


r 0 = q + 


dr 
dtj o 


J_ (—) 


2 \dt 2 jp 

dr\ 
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k 2 s 

At 2 = q + —. — At 2 
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d 2 r\ . 
— A t 
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At 2 + .= k*. — At — — — (3 s + 1) At 3 + . . . 

q2 " “ 


m =[-\ +(-) a t+ u 

{dt 2 )o [dfjp \dt 3 ) P 2 \ 

Das letzte Glied kann übrigens auch nur bei kleineren Periheldistanzen merklich werden. 
Mit diesem Wert bestimmt man 
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dt 4 )p q 2 2 q° 


*f(a- = -H d ^\ 
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und 
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Das in l f — — j sonst noch berücksichtigte Glied 4- ~ f l (a) kann für diese Ausgangsepoche 

(d*r\ 

wegen — =0 nur ein sehr kleiner Betrag sein, der kaum Anlaß geben wird, eine Neubestimmung von 

\dt 3 Jp 

l f[a -^-] vorzunehmen. 


III. 


Der Gang der Rechnung stellt sich also folgendermaßen: 

Man bestimmt zunächst den Abstand A t der nächsten Epoche vom Perihel und die Größen 


woraus 

r 0 = q + -^~ Af 2 , 

w 

und weiter die Ausgangswerte 


hH k* „ 0 , 

c — —, v — — (1 4- 3 s), 


(—) —cM + —cv M 3 , 
\dt) o 6 



\ 2 J «\dt) a 2 UWo 


“/(«) = 


1 


0 ) ‘ 


\2\dt\ 


0 



folgen. Aus den beiden letzteren erhält man unmittelbar ll f (a + w), womit nun das angegebene Verfahren 
begonnen wird nach dem Formelsystem 


log 


!cPr 

U t\ 


0)611/3 

log /(a + <ü) = logjp 4- 0*0362 co 2 v 
o) 2 d 2 r 


r = o) 2 n f 4* 


12 dt 2 


woraus wieder IX / (a 4- 2 co) und damit r für dieses nächste Argument erhalten wird, wobei nur zu 
bemerken ist, daß hier und bei den weiteren Fortsetzungen die Größe v bestimmt wird aus 


Es soll noch hervorgehoben werden, daß zur Bestimmung von F unter allen Umständen eine fünf¬ 
stellige Rechnung genügt. 

Dieses — wie man sieht, sehr einfache -— Verfahren der sukzessiven heliozentrischen Distanzen 
wird besonders im Falle großer Exzentrizitäten eine wesentliche Abkürzung der Rechnung bedeuten. 

Es soll dasselbe nun an einem Beispiel gezeigt werden. Ich wähle dafür die Erscheinung 1892 des 
periodischen Kometen Winnecke, dessen Exzentrizität von über 0*7 die Auflösung der Keppler’schen 
Gleichung schon recht unbequem werden läßt, andrerseits aber selbst in der Nähe des Perihels noch 
keinen Anschluß an parabolische Verhältnisse gestattet. 


5 * 
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Dem Teil der Ephemeride, der sich auf die Umgebung des Perihels bezieht, hat der damalige 
Berechner v. Haerdtl folgendes Elementensystem zu Grunde gelegt: 


woraus folgt: 


Epoche und Oskulation: 1892 Juli 4-0 m. Z. Berlin 
M— 0° 31' 14"73 


■k ~ 276° 11' 4 W 49 
ft = 104° 4' 37 w 05 
i= 14 D 31'33"64 
<p = 46° 33' 4' / 81 
(A =r 609"672 1260 


m. Äqu. 1890-0 


log a — 0-509 9402 
log p — 0-184 7438 


q — 0-886 5542 

Perihelzeit r=Juli 0-925007 m. Z. Berlin. 


Für die beiden Hilfsgrößen c und v 0 findet man 

löge = 6-43 6683, logv 0 = 7• 13 0197 

Soll die Ephemeride mit Juli 0 - 5m. Z. Berlin beginnen, so ist A t — —0-425007 und man erhält für 
diese Epoche 

r 0 = 0-886 5789 
| = — 0-000 1161 65 

o 


’dr 

dt 


—\ =+ 0-000 2732 94 

dt 2 Jo 

Daraus ergeben sich die Anfangswerte der beiden Summenreihen 

if[ a -^-\ = — 0-000 1947 0 


"/(*) - 


0-221 6219 0 


wobei für die Ephemeride ein zweitägiges Intervall angenommen wurde; daher auch 

0*0362 (o 2 v = 0*145 v. 


Mit diesen Werten kann nun die Rechnung begonnen werden, die sich für den Verlauf des Monates 
Juli nach dem angegebenen Schema so stellt, wie die beiden folgenden Tafeln zeigen. Es soll dazu nur 
nur noch bemerkt werden, daß dabei eine überflüssige Genauigkeit angewendet wurde; es würde eine 
Summierung auf sieben Stellen und eine vierstellige logarithmische Rechnung völlig ausreichen, um r 
auf sechs Stellen sicher zu erhalten. 

Ich möchte schließlich noch hervorheben, daß die einzelnen Posten dieses .Rechenschemas durch 
die einfachsten Operationen gewonnen werden, also immer unmittelbar hinzuschreiben sind, so daß die 
Raschheit der Ausführung desselben nichts zu wünschen übrig läßt. 
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(^Winnecke 1892 

2*5 

4*5 

6-5 

8'5 

Juli 

v. 

— 0*00135 

— 0*00137 

— 0*00I32 

— 0 * 00 I 3 I 

n /. 

0*221 7005 3 

0*222 0519 8 

0*222 6744 I 

0*223 5 6 45 1 

0 = tu 2 11/ . . . . 

0*886 8021 

o*888 2079 

0*890 6976 

0*894 2580 

1 d^r 

3 dß 

+ 909 

+ 903 

+ 892 

+ 877 

log c. 

9-94783 

9-94851 

9*94973 

9 * 95 I 46 

P-° . 

0-64337 

0*64196 

0-63947 

0 * 6359 I 

log(>—0). 

9*80846 

9-80751 

9*80582 

9*80340 

log Iß (p —0). 

6*27962 

6*27867 

6 * 27698 

6* 27456 

log o 3 . 

9'84349 

9-84553 

9-84919 

9*85438 

log F . 

6-43613 

6 ' 433 j 4 

6-42779 

6*42018 

0 * 145 v .... 

— 20 

— 20 

- 19 

— 19 

cßr 

dt 2 

o'ooo 2728 6 

o*ooo 2709 8 

o*ooo 2676 7 

0*000 2630 2 

r . 

0*886 893 

0*888 298 

0*890 787 

0*894 346 

(^Winnecke 

10*5 

12*5 

14*5 

1 

16*5 

Juli 

v. 

— 0*00130 

— 0*00127 

— 0*00123 

— 0*00118 

«/. 

0*224 7176 3 

0*226 1278 6 

0*227 7881 8 

0*229 6905 5 

O = (JU 2 II f . . . . 

0*898 8705 

0*904 5114 

0*911 1527 

0*918 7622 

1 d^r 

3 dß 

+ 857 

+ 833 

+ 807 

+ 777 

log 0 . 

9 ‘ 9537 o 

9*95641 

9*95959 

9*96320 

P-* . 

0*63130 

0*62566 

0*61902 

o*6i141 

log {jp 5 ). 

9*80024 

9*79634 

9*79170 

9*78634 

log Tß {p — 0) . 

6*27140 

6*26750 

6*26286 

6*25750 

log . 

9* 86110 

9*86923 

9*87877 

9*88960 

log F . 

6*41030 

6*39827 

6*38409 

6*36790 

0 * 1 45 v .... 

- 19 

- 18 

- 18 

- 17 

d*r 

dß 

0*000 2571 I 

o*ooo 2500 9 

o*ooo 2420 5 

0*000 2332 0 

r . 

0*898 956 

0*904 595 

0*911 233 

0*918 840 
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„ 

(j Winnecke 1892 

1 8 * 5 

20-5 

22-5 

24*5 

Juli 

v . 

— 0 * OOII2 

— 0-00108 

— 0*00102 

—0-00097 

“/. 

0*231 8261 2 

0*234 1853 4 

0*236 7581 I 

o-239 5339 6 

ci 

II 

c 

°'9 2 7 3045 

o*936 74H 

0*947 0324 

0-958 1358 

1 d*r 

+ 745 

+ 712 

+ 677 

+ 64I 

3 dt 2 





log 0. 

9-96722 

9-97162 

9-97637 

9*98143 

7 7 “ 0 . 

0*60287 

0*59343 

0-58314 

0*57204 

\ og ( p — 6 ) . 

9-78023 

9*77337 

9-76577 

9*75743 

log * 2 (jP —0). 

6-25139 

6*24453 

6-63693 

6-22859 

log G 3 . 

9•90166 

9 * 91486 

9-92911 

9*94429 

log F . 

6-34973 

6*32967 

6*30782 

6 * 28430 

0 * 145 v . . . . 

— l6 

- 16 

- 15 

- 14 

d 2 r 

o'ooo 2236 5 

o*ooo 2135 5 

0*000 2030 8 

o*ooo 1923 8 

dP 





r . 

0-927 379 

0*936 813 

0*947 100 

0*958 200 

$ Winnecke 

26-5 

28*5 

30*5 

i*5 

Juli 

August 

v . 

— 0 * 00090 

—o-ooo86 

— 0*00081 

— 0*00076 

«/. 

0*242 5021 9 

0-245 6519 7 

0*248 9724 9 

0*252 4530 3 

3 = u,2 11/ . . . . 

0*970 0088 

0-982 6079 

0*995 8900 

1*009 8121 

1 dir 

+ 605 

+ 569 

+ 533 

+ 498 

3 dP 





log g. 

9-98678 

999238 

9-99821 

0*00424 

P~ G . 

0*56016 

0-54756 

0-53428 

0*52036 

logO-o). 

9'7483i 

9-73843 

9-72777 

9*71630 

log 1P {jp 0). 

6*21947 

6-20959 

6-19893 

6*18746 

log G 3 . 

9*96034 

9-977H 

9'99463 

0*01272 

lOg-F. 

6-25913 

6-23245 

6 * 20430 

6*17474 

0*145 v . . . . 

- 13 

— 12 

— 12 

— 11 

d~r 

0*000 1815 5 

o-ooo 1707 4 

0000 1600 2 

0*000 1495 0 

dP 





r . 

0*970 069 

0*982 665 

o*995 943 

0*009 862 
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zu 

ß 

Juli 

d- V 

f~ä & 

V 

uz 






— o-ooo 1947 0 





°‘5 

o-ooo 2732 9 


0'22I 6219 4 



~ 43 



4- o-ooo 0785 9 



- 145 


2’5 

2728 6 


0*221 7005 3 



188 



0 35 H 5 



143 


4’5 

2709 8 


0*222 0519 8 



33 i 



0 6224 3 



134 


65 

2676 7 


0*222 6744 I 



465 



0 8901 0 



126 


8-5 

2630 2 


0-223 5645 I 



59 i 



1 I 53 1 2 



111 


IO 5 

2571 1 


0-224 7176 3 



702 



1 4102 3 



102 


IS S 

2500 9 


0*226 1278 6 



804 



1 6603 2 


(^Winnecke 

81 


i 4‘5 

2420 5 


0-227 7881 8 

1892 


1-0 

00 

00 



1 9023 7 



70 


16-5 

2332 0 


0-229 6905 5 



955 



2 1355 7 



55 


18*5 

2236 5 


0-231 8261 2 



1010 



2 3592 2 



37 


20*5 

2135 5 


0-234 1853 4 



1047 



2 5727 7 



23 


22 * 5 

2030 8 


0-236 7581 1 



1070 



2 7758 5 



- 13 


24-5 

1923 8 


o '239 5339 6 



1083 



2 9682 3 



4- 2 


26’ 5 

1815 5 


0-242 5021 9 



1081 



3 1497 8 



9 


28-5 

1707 4 


0-245 6 Si 9 7 



1072 



3 3205 2 



4- 20 


3 °’ 5 

1600 2 


0*248 9724 9 



— 1052 



4- 0-003 4805 4 





32*5 

0-000 1495 0 


0-252 4530 3 


Um nun ein Urteil über die Brauchbarkeit dieser Methode bilden zu können, seien hier die 
Logarithmen der so gefundenen Radienvektoren sechsteilig angegeben und nebenan die letzten zwei 
Stellen von log r, wie sie eine sechsteilige Rechnung nach der gewöhnlichen Methode ergibt. 


1892 log r 


Juli 0 - 5 

9-947717 

16 

2-5 

947871 

69 

4-5 

948559 

59 

6-5 

949774 

75 

8-5 

951506 

05 
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1892 

log r 


Juli 10 - 5 

9-953739 

38 

12-5 

956454 

56 

14-5 

959630 

32 

16-5 

963240 

41 

18-5 

967257 

57 

20-5 

971654 

55 

22-5 

976396 

98 

24-5 

981456 

56 

26-5 

986803 

02 

28-5 

992406 

06 

30-5 

9-998235 

36 

32-5 

0-004262 

62 


Die Übereinstimmung ist eine vollkommene, da die Abweichung der Größe und dem Sinne nach 
durchaus den Charakter der zufälligen Unsicherheiten einer sechsteiligen Rechnung tragen. 


IV. 

Aus r läßt sich nun die wahre Anomalie bestimmen; allerdings nicht direkt aus der Kegelschnitts¬ 
gleichung, da die Bestimmung von v aus r bei kleinen Exzentrizitäten überhaupt, bei größeren aber in 
der Umgebung des Perihels sehr unsicher wird. Es läßt sich aber aus dem Flächensatz 

dv j \\J p 

dt r 2 


und daraus v selbst mittels numerischer Integration mit völliger Schärfe bestimmen. 

Wjj 


Setzt man 


i/, so ist 


1 ra + (i -f i)' 


1 r a + 
Ja + i 


fdt=f 


a+\t + — 


12 


a + \t —— j oj 


ü -\r [i -] W 

_ÜL/n 

a + 


\ 2j J 

12 


l 2; J 


Da die beiden ersten Glieder hier vollkommen ausreichend sind, so erhält man als Differenz von 
zwei aufeinander folgenden wahren Anomalien » 

A vi — co f 

Was den Ausgangswert v 0 anbelangt, so findet man aus 

1 Z 1 2 ff 3 

— //dt = n f (a) + —f\a) + —/ n («) + ....+ Konst, 

w J 2 6 

wenn das Intervall zwischen der Ausgangsepoche und dem Perihel wieder mit A t bezeichnet wird, so daß 
A t 

n — — und Konst. = 0 ist: 
co 

-ffdt = -f{a) + ± . V!/i (a) + V . V 3 /“ («) + ••• 

co J co 2 co 2 6 co 3 













Ephemeridenrechnnng mittelst numerischer Integration . 


29 


demnach 

*0 = A tf{a) + d- A Pf (a) + -±- A Z 3 / 11 («)+.... 

2 (o 

Durch sukzessive Addition der obigen Größen Lv ergeben sich dann die wahren Anomalien für die 
folgenden Ephemeridenargumente. 

Nimmt man das obige Beispiel, die Erscheinung 1892 des periodischen Kometen Winnecke, wieder 

d v 

auf, so erhält man aus den gefundenen Radienvektoren zunächst folgende Werte für —, die gleich in 

d t 

Bogenmaß ausgedrückt sind: 


1892 

dv 

ß 

/ll 

/m 

dv 

Juli 

dt 




dt 



— 

— 

1 


0-5 

5584''03 

3 ? 99 



i° 33’ 4 r °3 

2-5 

5580-04 

17-65 

13*66 

o ’ ; 27 

1 33 0-04 

4* 5 

5562-39 

31 -04 

i3‘39 

0*48 

1 32 42*39 

6-5 

5531-35 

43-95 

12*91 

0*72 

1 32 11*35 

8-5 

5487-40 

56-14 

12 * 19 

0*85 

1 31 27*40 

IO * 5 

5431-26 

67-48 

n’34 

0*94 

1 30 31 * 26 

125 

5363-78 

77-88 

10*40 

1*12 

1 29 23*78 

14*5 

5285-90 

87-16 

9*28 

1 * 16 

1 28 5*90 

1 i6 ‘5 

5198-74 

95-28 

8-12 

1 -10 

1 26 38*74 

18*5 

5103-46 

I 02’30 

7*02 

1 * 29 

1 25 3*46 

20*5 

Ln 

O 

O 

O' 

108’03 

5*73 

1 • 26 

1 23 21 * 16 

22 ’ 5 

4893-13 

112*70 

4*67 

1 • 10 

1 21 33*13 

24^ 

4780-43 

Il 6’27 

3*57 

1 ■ 02 

1 19 40*43 

26*5 

4664-16 

Il 8*82 

2*55 

o*99 

1 17 44*16 

28*5 

4545-34 

I 20*38 

1*56 


1 15 45*34 

3°* 5 

4424-96 




1 13 44*96 
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Für den Ausgangswert v 0 erhält man, da A t = — 0-425007, 

A tf(a) = — 39' 33"25 

— \t 2 f l (a) = - 0"13 

2 co 

— A/ 3 / Il («) = - 0"04 

6 co 2 

demnach v 0 = — 0° 39' 33 // 42 

Die sukzessiven Differenzen Au und die daraus sich ergebenden v sind in nachstehender Tafel 
gemäß dem obigen Rechenschema zusammengestellt. 


1892 

Juli 

- 

0)/ 

1 

— ü>/II 

12 

V 

kv 

V 0 

°*5 

H- 3 ° 

6 1 4 v o 7 

- o< 

3 39 ’ 

33 r 42 

33 ? 4 



H- 2 • 28 

“ 4 - 3 

6 

6*35 


2-5 

3 

5 42*43 

2 

26 

32-93 

32*4 



-1- 2-25 

4 - 3 

5 

44*68 


4’5 

3 

4 53*94 

5 

32 

17*61 

17*4 



H- 2 * 19 

4 - 3 

4 

56-13 


65 

3 

3 33*75 

8 

37 

13-74 

13-2 



-4-2-09 

4 - 3 

3 

40*84 


8-5 

3 

1 58-66 

11 

40 

54*58 

54*2 



-h 1 ‘ 9 ^ 

4 - 3 

2 

0*62 


10-5 

2 

59 55*04 

14 

42 

55 ‘ 2 o 

54*6 



H- 1*81 

4 - 2 

59 

56-85 


12-5 

2 

57 29-68 

17 

42 

52-05 

5i*4 



+ 1-64 

4- 2 

57 

3 1 ’ 3 2 


14-5 

2 

54 44*64 

20 

40 

23-37 

22 * 1 



H- i *45 

4- 2 

54 

46*09 


i6-5 

2 

51 42-20 

23 

35 

9-46 

8-i 



-h 1 *26 

-4- 2 

5i 

43-46 


iS-5 

2 

48 24-62 

26 

26 

52-92 

51-8 



-4- 1 -06 

4- 2 

48 

25*68 


20*5 

2 

44 54*29 

29 

15 

i8*6o 

17 • 2 



-h 0*87 

4- 2 

44 

55-16 


22 • 5 

2 

4 i I3*56 

32 

0 

13-76 

12*6 



+ 0-69 

4- 2 

41 

14-25 


24-5 

2 

37 2 4*59 

34 

41 

28-01 

27 • 6 



-4-0-51 

4- 2 

37 

25 • 10 


26*5 

2 

33 29-50 

37 

iS 

53*ii 

5 2 *9 



H-o-34 

4- 2 

33 

29*84 


28*5 

2 

29 3 °* 3 ° 

39 

52 

22-95 

22 • 1 



+ o-18 

4- 2 

29 

30-48 


30-5 



42 

21 

53-43 

52-1 












Epheineridenrechnung mittelst numerischer Integration . 


31 


In der letzten Kolumne sind die Sekunden von v angegeben, welche bei einer sechsstelligen 
Rechnung nach der gewöhnlichen Methode mittelst der exzentrischen Anomalie erhalten wurden. 

Es zeigen sich hier nicht unbeträchtliche Abweichungen in den nach beiden Methoden resultieren¬ 
den Werten für die wahre Anomalie. Dieselben rühren aber zum größten Teil von der Unsicherheit der 
sechsstelligen Rechnung nach der gewöhnlichen Methode her. DieErmittlung der wahren Anomalie aus 
der mittleren bedingt bei großen Exzentrizitäten eine gewisse Unsicherheit, die schon darin ihren Aus¬ 
druck findet, daß man aus kleinen Werten M relativ große Beträge v zu ermitteln hat. 


Aus 


folgt in der Tat 


r 2 dv — !z\J p dt 


dv— — \J 1-s 2 dMz 


(1 + £ COS V ) 2 
~—(\— s 2yi z ~ 


dM. 


Für den vorliegenden Wert der Exzentrizität ist der Faktor von dM im Perihel ungefähr 9*2 und 
wird erst bei v = 110° der Einheit gleich. In der Umgebung des Perihels können daher Abrundungsfehler 
in M sehr merklich in v eingehen. 

Die hier vorgeschlagene Methode der numerischen Integration, die nicht aus dem Fortschreiten 
dv 

von M, sondern direkt aus — das Fortschreiten der wahren Anomalie bestimmt, ist naturgemäß von 

dt 

dieser Unsicherheit frei. Es kann dies aus irgend einem der obigen v ersehen werden. Prüft man darauf¬ 
hin etwa die beiden letzten ziemlich stark voneinander abweichenden v für Juli 30*5 

42° 21 ; 53"4 aus der numerischen Integration und 

52 // l aus der gewöhnlichen sechsstelligen Rechnung. 

Die entsprechende mittlere Anomalie ist 

M— 5° 0' 31 "04. 

Die gewöhnliche sechsstellige Rechnung, für die Af — 5° O'SF'O anzunehmen ist, ergibt für v den 
obigen zweiten Wert; aus einer siebenstelligen Rechnung, die noch Hundertelsekunden in M berück¬ 
sichtigen kann, resultiert hingegen 

v — 42° 21'53"1, 

ein Wert, dessen Abweichung vom vorigen der Größenordnung nach mit der obigen Überlegung über¬ 
einstimmt und der dem aus der numerischen Integration gefundenen beträchtlich näher liegt. Die noch 
übrig bleibende Differenz hält sich innerhalb der naturgemäßen Unsicherheit einer längeren Serie von 
numerischen Integrationen. 

Da demnach bei erhöhter Genauigkeit der Aufwand an rechnerischer Arbeit beträchtlich kleiner ist 
als bei einer sechsstelligen Rechnung nach der gewöhnlichen Methode, so ist auch bei der Ermittlung 
von v — wenigstens im Fall großer Exzentrizitäten — der Vorzug dieser Integrationsmethode wohl 
unzweifelhaft. 


5* 
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V. 


Man wird nach Ermittlung der Größen r und v in der gewöhnlichen Weise mitHilfe der Gauß’schen 
Konstanten a , b, c, A ', B', Ö die rechtwinkligen heliozentrischen Koordinaten rechnen. Es liegt allerdings 
der Gedanke nahe, nachdem r bestimmt ist, diese aus den Differentialgleichungen nach ähnlichen Methoden 
direkt zu erhalten. Man kommt tatsächlich wieder auf sehr einfache Relationen; trotzdem ist es mit Rück¬ 
sicht auf die Raschheit und Bequemlichkeit des gewöhnlichen Verfahrens sehr fraglich, ob durch die Aus¬ 
dehnung der angegebenen Methode auch auf diesen Teil der Ephemeridenrechnung ein ökonomischer 
Vorteil erzielt wird. 

Die Zusammenstellung der dynamischen Gleichung 

^f + i>A = 0 

dt 2 r 3 


mit der Formel der numerischen Integration 


1 ur 1 d * x 
— x — u f q- 


1 


O) ' 


12 dt 2 240 




ergibt 


1 f n 


X—(ö 2U f 


240 u f 


1 q-o) . — 


12 r* 

wofür man auch, wenn nicht außergewöhnliche Verhältnisse (sehr kleine Periheldistanzen oder sehr große 
Intervalle <*>) vorliegen, setzen kann 

k 2 O) 2 

r3 


X = <ü 2U f tl — 


12 


d 2 x 


so daß sich aus der für — 1 anzulegenden Summentafel die rechtwinklige Koordinate direkt und in sehr 

dt 2 

einfacher Weise rechnen läßt. 

Die Ausgangswerte für die beiden Summenreihen lassen sich leicht aus den Gauß’schen Konstanten 
finden. Es ist 

x 0 = r Q sin a sin {A! q- v 0 ) 

= ksm J* [cos (A! q- v 0 ) q- s cos A '] 
dt)o \/p 


und 


’d 2 x\ _ k 2 x 0 


2) iü\dt) 0 2\dt*) 0 


12 




wodurch die Ausgangswerte 
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und 


ü m=\x o 


12 ydtjQ 


1 

240 


-/ n («) 


gegeben sind. 

Die Relationen sind, wie man sieht, zwar auch sehr einfach, doch kaum derartig, daß man Anlaß 
Hätte, von der gewöhnlichen Weiterführung der Ephemeridenrechnung abzugehen. 

Der Schwerpunkt der hier auseinandergesetzten Methode bleibt in der gleich bequemen Ermittlung 
von r und v, wie groß auch immer die Exzentrizität sein mag, ein Verfahren, das sich insbesonders 
bei einigermaßen beträchtlichen Werten von s der gewöhnlichen Methode gegenüber sehr im Vorteil 
befindet. 

Ich möchte zum Schluß auf den nicht unwichtigen Umstand hinweisen, daß die vorgeschlagene 
Methode eine sehr leichte durchgreifende Kontrolle gestattet: Wenn man nach einer Serie von Werten 
r und v für den letzten Ort die gewöhnliche direkte Methode in Anwendung bringt, so verifiziert ein 
übereinstimmendes Resultat die ganze bis dahin geführte Rechnung. 
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